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ABSTRACT 

AK Sco is a unique source: a ~10 Myrs old pre-main sequence spectroscopic 
binary composed of two nearly equal F5 stars that at periastron are separated 
by barely eleven stellar radii so, the stellar magnetospheres fill the Roche lobe 
at periastron. The orbit is not yet circularized (e=0.47) and very strong tides 
are expected. This makes of AK Sco, the ideal laboratory to study the effect 
of gravitational tides in the stellar magnetic field building up during pre-main 
sequence (PMS) evolution. In this letter, the detection of a highly disturbed 
(cr ~ 100 km s~^) and very dense atmosphere (ue = 1.6 x 10^° cm~^) is reported. 
Significant line broadening blurs any signs of ion belts or bow shocks in the 
spectrum of the atmospheric plasma. The radiative loses cannot be accounted 
solely by the dissipation of energy from the tidal wave propagating in the stellar 
atmosphere; neither by the accreting material. The release of internal energy 
from the star seems to be the most likely source of the plasma heating. This is 
the first clear indication of a highly disturbed atmosphere surrounding a pre-main 
sequence close binary. 

Subject headings: stars: pre-main sequence, stars: magnetic fields, binaries: spec- 
troscopic 



1. Introduction 

Pre-main sequence (PMS) stars are surrounded by very powerful magnetospheres that 
extend 2-4 stellar radii to the limits of the inner disk boundary. The magnetospheres are 
powered by the stellar magnetic field, the accreting infalling material and the shearing at 
the boundary between the magnetosphere and the inner disk (see Gomez de Castro 2009 
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for a recent review). Magnetospheric heating processes during pre- main sequence evolution 
are poorly studied, as well as the relative relevance of the contributions mentioned above. 
This is partially caused by the similar densities (10^ — 10^^ cm~^) and temperatures (lO^'^ — 
10^-^K) of several distinct physical components in the circumstellar environment (the disk- 
star shear layer, the outflow and the accretion flow). Moreover, the high temperature of 
the plasma makes thermal broadening very large (~ 70km s~^) leaving ultraviolet high 
spectral resolution monitoring as the only technique to resolve the various contributions 
since variability time scales are expected to be very different for the various components (see 
i.e. Gomez de Castro & Verdugo 2007). 

Close pre-main sequence binaries are ideal targets to unveil some of the magnetospheric 
properties. Circumstellar accretion disks are difficult to accommodate in these systems 
since the strong gravitational tides break them into rings. The inner boundary of the cir- 
cumbinary disk is pushed outwards by the binary orbit thus the magnetosphere of each 
individual star is not locked to the disk rather, it is expected that a strong coupling between 
the gravitational tides and the stellar magnetic fields ought to be a major source of atmo- 
spheric/magnetospheric heating. The investigation of the atmospheric properties of close 
pre-main sequence stars provides important clues regarding the possible role of magneto- 
spheric dynamos in pre-main sequence evolution. 

AK Sco is a unique target for this purpose; it is the only pre-main sequence binary 
known to date with small eccentricity and composed of two nearly equal stars that get as 
close as 11 R* at periastron passage. No signs of enhanced accretion have been detected 
at periastron (Andersen et al. 1989, Alencar et al. 2003). Polarimetric observations have 
reported the largest variations of percent polarization and position angle detected in a pre- 
main sequence binary (Manset et al. 2005). Variations are found to be periodical with 
period equal to the gravitational tide period (half the orbital period). The main parameters 
of the binary are summarized in Table [H In this letter, the profiles of the (semiforbidden) 
intercombination lines of the optically thin magnetospheric tracers: the CIII] at 1908A 
and the Si III] at 1892A spectral lines, are analyzed. From the analysis, it is concluded 
that there is an extended region of radius ~ 1.3 — 4.9R* around each star containing hot 
(logTg = 4.8) and highly perturbed gas that produces a gaussian-like profile broadening with 
cr ~ 100 km s~^. Although this broadening is highly suprathermal it could be sufalvenic if 
AK Sco's magnetic field is similar or higher than that reported for other pre-main sequence 
stars (see i.e. Johns-KruU, 2007). 
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2. Observations and measurements 

AK Sco was observed with ttie Space Telescope Imaging Spectrograph (STIS) on March 
10, 2001. The observations were carried out with grating G230M to target the semiforbidden 
lines of Si III] at 1892A and C III] at 1908A. The spectra were obtained at phase 0.2814 
according to Andersen et al (1989) ephemeris (see Table 1). At this phase, the radial velocity 
of both components is the same thus the profiles represent the co-added contribution of both 
components in the same velocity space. As all indications are that the two components are 
equal in all respects, including spectral type, the contribution from both stars is expected 
to be the same. 

AK Sco displays strong and very broad Si III] and C III] profiles. The observed line 
profiles have been converted to velocity distribution using the Doppler shift formula. Data 
have been processed with the Routine Science Data Pipeline (RSDP). The major source of 
inaccuracy in the calibration is the centering of the target in the aperture that can account 
for as much as 2.8 km s~^. The exposure duration was 1944 s that were splitted into three 
subexposures; the average value for each pixel and the deviation from this average are plotted 
in Fig. 1. A curve fit to a single gaussian plus a constant is also shown. The dispersion of the 
gaussians are a{C III])=109.5 km/s and a{Si III]) = 106.8 km/s. The value of x^/^ is 3.6 and 
1.3 for this best gaussian fit to the Si III] and C III] lines, respectively. The instrumental 
profile of HST/STIS with the grating G230M has a FWHM of 30 km/s and the stellar 
rotation velocity is vsini = 18.5 km s^^, thus line broadening is more than one order of 
magnitude larger than expected from the rotational broadening of the atmospheric lines. 

The lines fluxes are: F(Si III]) = (4.2 ± 0.1) x lO^^^g^g/s/cm^ and F(C III]) = (2.9 ± 
0.1) X 10~^^erg/s/cm^, after correction for extinction using Valencic et al. 2004 UV law (R= 
Av/EB-y=4.3 from Manset et al. 2005 and Ay = 0.5 ± 0.1 mag from Alencar et al. 2003). 
The lines ratio is calculated to be F(Si III])/F(C III]) = 1.4 ± 0.2. This value is similar 
to that derived from earlier (1982 and 1988) low resolution observations obtained with the 
International Ultraviolet Explorer (lUE) (Gomez de Castro & Franqueira 1997). From the 
lUE low dispersion spectra a C IV(AA1548, 1550A) flux of (1.6 ±0.1) x 10~^^erg/s/cm^ and 
a ratio F(C IV)/F(Si III]) of 1.9 ± 0.4 was derived. There are only three low resolution 
observations carried out with the lUE, two of them during the same cycle at periastron 
(phase=0.123) and at apastron (phase=0.708). No significant variations in the C III], Si III] 
and C IV fluxed are found though the lUE fluxes are about a factor of ~ 2 higher than that 



iThe FWHM of the profiles is FWHM(Si III])=251 km/s and FWHM(C III])==258 km/s. 

^There are liowever, variations by a factor of ~ 2 in the O I (A1304A) and C II resonance multiplet at 
I335A; these lines trace cooler and lower density plasma. 
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measured with the HST. 



3. 



Constraints from the lines flux 



Assuming coUisional equihbrium, the F(Si III])/F(C III]) and F(C IV)/F(Si III]) ra- 
tios indicate that the physical conditions in this Hot Line Emission Region (HLER) are 
logne(cm~^) = 10.2 and logTe(K) = 4.8. An optically thin coUisional plasma with these 
properties radiates 2.4 x 10~^ of the total flux in the C III] line; thus the luminosity of 
the HLER is 3.1xl0'^^erg/s or O.O8L0 (0. 041^0 per star assuming the same contribution per 
component); a distance, d, of 145 pc to AK Sco has been assumed. To account for this high 
luminosity an accretion rat^ of Ma = 0.9 x 10~^ MQyr~^ is required. The viscous release 
of this gravitational energy in the circumbinary disk would affect signiflcantly the infrared 
spectral energy distribution (SED) however, good flts are obtained with passive disk models 
suggesting that the accretion rate is not very much above 1 x 10"^ M0yr~^ (see i.e. Alencar 
et al. 2003 for a discussion on the SED flt). 

The characteristic size of the HLER, Re, is given by. 



where Rg is the radius of an equivalent homogeneous spherical envelope around each star, 
with < ji > representing the average emissivity in the (Si HI]) lines formation region, which 
for this study, has been assumed to be of the order of ji = 2.38 x 10~^erg s~^ cm~'^ using 
the atomic parameter, from the Chianti Atomic Data B^g) . Here / represents the volume 
flUing factor of the radiating plasma. For / ~ 0.1, Re reaches 1.7 R*. This is an indication of 
the HLER being signiflcantly more extended than the atmosphere of a cool main sequence 
star. 



Line broadening is highly suprathermal (thermal broadening at 10^'^K generates a 
FWHM of 73 km/s, a factor of ~ 3 smaller than the observed). In principle, such a broad ki- 



■^A radius of 0.43 AU (three times the seminiajor axis of the orbit) is estimated for the gap of the 
circumbinary accretion disk based on numerical simulations of the evolution of binary systems (Artymovicz 
& Lubow, 1997). 

'^URLtwww. damtp.cam.ac.uk/user/astro/chianti/ 




4. Constraints from the line broadening 
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netic velocity distribution of the emitting ions could be produced, as in the Sun, by resonances 
between Alfven waves propagating in the atmosphere and the giro-frequency of the ions, driv- 
ing the Si III] and C III] ions to the very high kinetic temperature of Tkin = 0.75 x lO^K (sec 
i.e. Cranmer 2002). However, the strength of this broad component must correlate with the 
coronal X-ray surface flux produced by the bremsstrahlung radiation of the high velocity 
electrons, as shown for magnetically active stars by Wood et al. 1997. Though AK Sco 
X-ray surface flux. Fx — 6.49 x 10^ erg/s/cm^ (derived from Steltzer et al. 2006 based on 
CHANDRA observations) is similar to the detected in class IV-V F5 stars (see Ayres et al. 
1995), the CIV flux exceeds by ~ 4.5 orders of magnitude the flux predicted by the active 
stars CIV-X ray surface flux correlation. Thus, the HLER radiation cannot be accounted 
for solely by this mechanism. 

The high symmetry of the profiles is suggestive of them being formed in a rotating 
ionized ring or belt such as the observed around RW Aur (Gomez de Castro & Verdugo 
2003). However the characteristics of AK Sco orbit imposes severe constraints on the possible 
location of the material with Keplerian-like orbits. On the one hand, circumstellar material 
is forced to be between the stellar corona (~ 2i?^,) and the location of the Lagrange point 
at periastron (5.5i?*). Plasma in such a hypothetical ring systems surrounding each star 
would be orbiting with velocities in the range 180 < King < 280km s~^ and thus, could 
only contribute to the line wings. On the other hand, the binary cleans the inner region of 
the circumbinary disk to about 0.43 AU thus, any ionized material in this boundary would 
be rotating at Keplerian velocities ~ 75km s~^. As a consequence, the profile broadening 
can neither be accounted for by rotating rings around the stars nor by a circumbinary 
belt. Moreover the low rotation velocity of the stars precludes that line broadening can 
be accounted by an extended magnctosphere, rigid-body like rotation as shown in Fig 2. 
Line profiles have been simulated for various emissivity and density profiles for a corotating 
magnctosphcric shell and belt. The energy radiated in the lines per unit volume, e^ad,;, at a 
given distance, r, from the center of the star is given to first order as: e^ad,; = Ht) x p{r) oc 
where j(r) represents the line emissivity and p the density of emitting ions. As shown in Fig 2, 
simple spherical and cylindrical distributions (simulating a shell-like or a belt-like structure, 
respectively) cannot account for the observed widths of the HLER fines. The optimal fit is 
obtained by adding a Maxwellian velocity distribution with variance from a —74 km/s for 
a corotating ring peaking at d=5.5i?* (x^ = 27.5 and x^/^ = 1-44) to a =88 km/s for a 
spherical envelope with e^ad oc {x^ = 19.5 and x^/^ = 1.2). In all cases, the Maxwellian 
velocity field is highly supersonic and blurs the possible kinematical signatures of rotating 
structures around the star. 
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5. Conclusions: on the source of the HLER heating 

AK Sco has an UV excess that accounts by as much as 0.04 Lq per component. This 
energy cannot be fed into the system solely by the gravitational energy lost by the accretion 
flow channelled by the surrounding disk since the accretion rates required are about one 
order of magnitude larger than that derived from the infrared properties of the disk. Even if 
the contribution of the fresh deuterium to the nuclear reactions is taken into account (Siess 
et al 1997), the accretion rate is a factor of ten smaller than that required to produce the 
observed luminosit}0. 

The gravitational energy stored in the epicyclic motion of the relative orbit is eo = 
i M-j^+M2 ^^'^^ ~ 2-36 X 10'^^ erg where Mi and M2 are the masses of the two components 
of the system, e the eccentricity, a the semimajor axis and uj = ^ with P the period. If 
the gravitational energy were dissipated through the HLER radiation, it would be damped 
in a time scale of: Tude = Co/O.OSLq = 2.4 x lO^yr"^ which is rather close to the age of 
AK Sco derived from the evolutionary tracks and the system should be already circularized. 
Henceforth, an additional source of energy needs to be hypothesized to power the HLER. This 
source also, should be acting in other pre-main sequence stars and be partially responsible 
of the observed UV excess. 

Finally, notice that the magnetic diffusivity in the HLER, rj, is very small (just 0.0045 
cm^/s) leading to very high Reynolds magnetic numbers and magnetic flux freezing. This has 
two important consequences: (1) the tidal deformation wave induces magnetic multipoles into 
the stellar dipolar flelds that can be responsible of the unusually strong periodic polarimetric 
variations of the system and (2) the Alfven waves induced by the tide will have very large 
damping lengths due to magnetic resistivity effects. The non-linear wave-wave interaction 
will lead to a turbulent cascade that will end at a scale determined by the dominant damping 
mechanism. For the HLER physical properties, this will be collisions between the ionized 
HLER material and the neutral infalling gas from the accretion flow that, in turn, strongly 
depends on the ill determined relative densities of these two components. If the observed 
broadening of the HLER lines is interpreted as a sign of magnetohydrodynamical turbulence 
this would imply that the inferred mean-fleld strength of the turbulent component is 3.5 — 
4.6G. This value is reasonable for the magnetosphere of a T Tauri star where typical surface 



^An estimate of the luminosity excess caused by accreting fresh deuterium can be made from the 
nuclear energy released per reaction is Qo = 5.5 MeV. The luminosity released is given by: L]j = 
{l/2)(rnD/'m)nAMaQD where rriD/m is the mass fraction of deuterium in the accreting matter, ua the Avo- 
gadro number and the factor of (1/2) is introduced to account for the Deuterium mass number. Assuming a 
mass fraction deuterium to hydrogen of 10^^, solar abundances and an accretion rate of 1.0 x 10^^ MQyr^^, 
the luminosity excess will be 0.004 Lq. 
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fields are 0.3-1 kG (Johns-KruU, 2007). However, this value would imply that the /3 of the 
atmospheric plasma is as high as 0.14!. 

To conclude, the observations presented in this letter indicate that AK Sco is the first 
pre-main sequence binary for which evidence of a highly perturbed atmosphere has been 
reported. However, conclusive evidence on the source of the perturbation and the mechanism 
that drives it, cannot be inferred solely from these data. Further UV observations and 
monitorings are required. 

This work has been partly financed by the Ministry of Education of Spain through 
grant AYA2007-67726 and by the Comunidad Autonoma de Madrid by grant CAM-S- 
0505/ESP/0237. 
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Fig. 1. — C III] (top) and Si III] (bottom) profiles of AK Sco obtained with the Hubble 
Space Telescope and the Space Telescope Imaging Spectrograph. Data have been processed 
with the Routine Science Data Pipeline (RSDP). The major source of inaccuracy in the 
calibration is the centering of the target in the aperture that can account for as much as 
2.8 km s~^. The exposure duration was 1944 s that were splitted into three subexposures; 
the average value for each pixel and the deviation from this average (error-bar) are plotted. 
A curve fit to a single gaussian plus a constant is also shown. 



- 10 - 





Fig. 2. — Predicted C III] line profiles for equatorial and spherical distributions of the 
emitting plasma. Top and bottom panels represent synthetic profiles for spherical-like and 
disk-like mass distributions. Three types of dependences of the line emissivity with the radius 
have been computed for the two possible mass distributions: erad,i oc with a = —2, and 
2, as indicated in the figures. The synthetic profiles are represented in the left panels. It is 
necessary to add a macro turbulence field to fit the profiles, as shown in the right panels. 
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Table 1: Main properties of AK Sco spectroscopic binary system'^ . 

Property Value 

Projected semimajor axis'' asini — 30.77 ± O.I2R0 

Eccentricity'' ^ e= 0.47 

Orbital period'' ^ P=13.609 d 

Inclination^ i = 65° - 70° 

Age= 10-30 Myrs 

Spectral typc^ <^ F5 

Stellar Mass^ M, = 1.35 ± O.OTM© 
Radius'^ = 1.59 ± 0.35Rq 

Projected rotation velocity'' '^^ vsini — 18.5 ± 1.0 km 



"The ephemeris is T{min) = F J£)2446666.380 + 13.609 x E from Andersen et al. 1989 
''from Andersen et al. 1989 
'^from Alencar et al. 2003 

"^The system is composed by two nearly identical stars so these parameters are the same for both components. 



